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(54) Method and device for digHally demodulating a frequency modulated signal 



(57) The present invention represents a new 
approach in the field of digitally demodulating a fre- 
quency modulated signal as it is used for example for 
fm-audio broadcasting. The inventive algorithm 
improves the signal to noise ratio of the demodulated 
signal for carrier to noise ratios below the fm-threshold. 

According to the present invention the spikes in the 
demodulated signal which are typical to occur for carrier 



to noise ratios below the fm-threshold are detected and 
eliminated. The detection is done by interpolation and 
other aiteria of the unwrapped phase tpu of the fre- 
quency modulated signal. The unwrapped phase Vu is 
calculated as the absolute phase which is not limited to 
the range [-n, ...,7i]. 
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Description 

[0001] The present invention is concerned with a method and a device for digitally demodulating a frequency modu- 
lated signal, for example a stereo broadcast signal, in particular with the extension of the fm-threshold of a digital fm- 
receiver. 

[0002] The modulating signal s(t) is frequency modulated: 



10 



afmW = COs|^CDct + AcD J s(T)dT j. 



(0) 



wherein coq denotes the carrier frequency and Aco denotes the frequency deviation. 

[0003] Figure 10 shows a state of the art digital fm-demodulator, which makes use of the CORDIC-algorlthm. The 
15 received fm-modulated signal af^W is down-converted to the IF (intermediate frequency) by an analog frontend block 
1 to an IF-signal b(t). This IF-signal b(t) gets sampled by an A/D-converter 2 with the sampling period T. 

b(kT) = B • cos CO ,pkT + Aca f s(T)dT , (1) 



20 



and converted to a complex baseband signal 



25 KT 

c-r(kT) = B •e^^"js(T)dT. (2) 



30 by an IQ-generation means 3. There are several algorithms existing for a digital fm-demodulation. The fm-demodulator 
shown In figure 1 1 makes use of the CORDIC-algorithm. The fm-demodulation is principally done by differentiating the 
angle (Pc(kT) of the complex baselsand signal Cj{WT): 



35 



40 



45 



kT 

cpc(kT) = Aco J s(T)dx. (3) 



= A© s(kT) 

kT 



[0004] With the approximation 



<Pc(kT) - q>c((k-l)T) 
T 

kT 



so follows: 



s(kT) « ^((Pc(kT) - <Pc((k-1)T)). (4) 

55 

wherefrom it can be seen that the fm-demodulation performed by the CORDIC block 4 and the differentiation block 7 is 
done by a determination of the angle of two samples of the complex baseband signal c-^kT) in the CORDIC block 4 and 
a subtraction in the differentiation block 7. 
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[0005] The CORDIC algorithm calculates the modulo 2n divided phase <p(kT) of the complex baseband signal c^kJ) 
by a transfbmfiation from cartesian coordinates to polar coorcfinates. 

(p(kT) = mod2^{<Pc(k"0} = modg^jA© J sCxjdxj (5) 



[0006] TWs CORDIC-algorithm is well-known. Rgure 1 1 and the following table 1 give a short introduction to this algo- 
10 rithm. Figure 1 1 shows the complex plane with the pointers of 5 samples of the complex baseband signal c-KkT). The 
argument range of the CORDIC-algorithm is limited to [-ti. .... n]. Th^-efore, the phase <p(kT) calculated with the 
CORDIC-algorithm is modulo 2n divided. The following table 1 shows the absolute phase (Pc(kT) of the complex base- 
band signal CjCkT) and the modulo 2n divided phase q)(kT) calculated with the CORDIC- algorithm. 



Table 1 



k 


1 


2 


3 


4 


5 


9(kT) 


2nJ5 


471/5 


-4n/5 


-7l/5 


7l/5 


<Pc(kT) 


2k/5 


471/5 


SrUS 


9n/S 


1W5 



[0007] State of the art demodulators using the CORDIC-algorithm demodulate the fm-modulated baseband signal 
C-KkT) by the subtraction explained in equation (4). Since the CORDIC-algorithm calculates the modulo 2n divided 
25 phase <p(kT) as explained above, the subtraction is canied out therewith and the following modulating signal s(kT) is 
obtained: 

s(kT) « ^(<P(k'n - q>((k-1)T)). (6) 

30 

[0008] The following table 2 shows the phase difference calculated within this operation for the modulo 2n divided 
phase (p(kT) and for the absolute phase 9c(kT): 

35 

Table 2 



k 


2 


3 


4 


5 


<p(kT)-9((k-1)T) 


271/5 


-Sti/S = -271 + 271/5 


371/5 


271/5 


<Pc(kT)-9c((k-1)T) 


2rd5 


271/5 


3n/5 


271/5 



[0009] It can be seen tfiat at the sample k = 3 for the phase difference of the nxxiulo 27c divided phase (p(kT) calculated 
with the CORDIC-algorithm an en^or of -27i occurs. In state of the art demodulators this error gets eliminated by a mod- 

45 ulo 2 71 operation carried out by an overflow in a fixed point implementation. Therefore, the valid range of the phase dif- 
ference qKkT) - <p((k-1) T) is limited in a fixed point implementation to [-tc. .... nj. So the value -Sn/S = -Zn + 2rc/5 causes 
an overflow and therefore the wrong value -Sn/S is modulo 27c divided to the correct value 2nJ5. 
[0010] The object underlying the present Invention is to improve the signal to noise ratio of the demodulated signal 
for canrier to noise ratios below the fm-threshold. 

so [001 1 ] According to the present invention, this object is solved by unwrapping the modulo 27i divided phase <p(kT) cal- 
culated during the transformation from cartesian coordinates to polar coordinates and the detection and elimination of 
phase jumps in the unwrapped phase to generate a corrected phase that is then used in the differentiation stage. 
[001 2] The method according to the invention is defined in independent claim 1 . Dependent claims 2 to 1 9 show pre- 
ferred embodiments of the inventive method. Devices that are capable of carrying out the inventive method are defined 

55 in claims 20 and 21 . 

[001 3] Further objects, features and advantages of the present invention will be better understood from the following 
analyzation of the state of the art demodulator and detailed description of preferred embodiments according to the 
present invention which are only to be seen as examples taken in conjunction with the accorrpanying drawing^. 
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wherein: 

Figure 1 shows the principal block diagram of an fm-demodulator according to an embodiment of the present 
invention; 

Hgure 2 shows the absolute phase cpcC^'O. i> e- the ideal phase having an infinite earner to noise ratio, and the 
phase cpu(kT) calculated by the transformation from the cartesian coordinates into polar coordinates in an 
unwrapped state with a carrier to noise ratio of 10 dB of a frequency modulated signal, here, the modulat- 
ing signal is a stereo-multiplex-signal with music information as it is used in fm-broadcasting; 

Figure 3 shows another unwrapped phase signal <pJ}(T) with a carrier to noise ratio of 10 dB of a frequency mod- 
ulated signal in an enlarged scale, here, the OKXIulating signal s(t) is again a stereo-multiplex-signal with 
music information; 

Figure 4 shows a flow chart for the generation of the unwrapped phase 9u(*^Tli 

Figure 5 shows the flow chart of the principal idea of the detection and elimination of the destortion of the 
unwrapped phase <Pu(I^T) shown in figure 2 and 3; 

Rgure 6 shows the unwrapped phase <Pu(kT) with a canrier to noise ratio of 10 dB and the phase q>d(kT) corrected 
according to the present invention; 

Figure 7 shows a flowchart of the blocks for unwrapping the phase, detection and elimination of the distortions and 
differentiation according to a second embodiment of the present invention; 

Figure 8 shows a direct realizatbn of the blocks for unwrapping the phase, detection and elimination of the distor- 
tions and differentiation according to the second embodiment of the invention; 

Figure 9 shows the signal to noise ratio at the demodulator output for a state of the art fm-demodulator and a fm- 
demodulator according to the present invention; 

Figure 10 shows a block diagram of a state of the art fm-demodulator which makes use of the CORDIC-atgorithm; 
and 

Rgure 11 shows the complex plane with the pointers of 5 samples of the complex baseband signal CT(kT). 

[001 4] Figure 1 shows the principal block diagram of a device for digitally demodulating a frequency modulated signal 
afm(t) according to an embodiment of the present invention. Blocks having the same functionality as those shown and 
described in connection with the prior art device shown in figure 1 0 are denoted with the same referencesigns. The dif- 
ferentiation block 7 of the state of the art device produces an overflow which is not allowed according to the present 
invention. Therefore, the device according to the present invention includes a differentiation block 6 not producing an 
overflow. Additionally to the prior art device the demodulator according to the present im^entbn comprises a digital 
threshold extension unit 5 receiving the phase (p(kT) from the CORDIC block 4 that calculates a corrected phase (p^ikT) 
which is then input to the differentiation block 6. An additional digital IF-filter can be placed in front of the CORDIC block 
4. The IF-filter is necessary for the neighbour channel suppression. 

[001 5] According to the embodiment of the present invention the digital threshold extension unit 5 comprises a phase 
unwrap unit 5a to calculate the unwrapped phase (Pu(kT) from the phase <p(kT) output of the CORDIC block 4 and a cor- 
rection unit 5b to detect and eliminate 2k phase jumps in the unwrapped phase (Pu(kT) and to output a corrected phase 
(|)e(kT) which Is then input into the differentiation block 6. The unwrapped phase <Pu(kT) Is calculated as the absolute 
phase which is not limited to the range [-n n]. 

[0016] The fm-demodulation done with the unwrapped phase <Pu(kT) without the correction unit 5b results in exactly 
the same values as the state of the art fm-demodulator described before. Ihis behavior is independent of the carrier to 
noise ratio (CNR) of the fm-signal. 

[0017] Figure 2 shows the absolute phase <PcQ^T) of a frequency modulated stereo-multiplex signal having a falling 
slope which is calculated with an infinite carrier to noise ratio. This absolute phase (Pc(kT) is shown with a dotted line. 
Figure 2 shows also the corresponding unwrapped phase 9u(kT) with a carrier to noise ratio of 1 0 dB. This unwrapped 
phase <Pu(kT) is shown with a solid line and it can clearly be seen that it includes a phase jump of 2n in regard to the 
absolute phase that is beginning with the value of k = 0. 

[001 8] The 271 phase jump of the unwrapped phase signal (Pu(kT) at sample k = 0 is typical for carrier to noise ratios 
below the fm-threshold. After the differentiation this phase jump produces a noise peak in the demodulated signal. Such 
a noise peak would also occur when the differentiation of the modulo 2n divided phase qp(kT) of the state of the art fm- 
demodulator is performed. These peaks in the frequency denruxlulated signal s(kT) are the reason for the dramatically 
decreased signal to noise ratio below the fm-threshold. 

[0019] These +2rt phase jumps of the unwrapped phase cPu(kT) appear in the normal mode of operation of the fm- 
receiver only at the falling slope of qpu(kT). Phase jumps with the value -2n appear only at the rising slope. 
[0020] According to the present invention such phase jumps are detected and eliminated to obtain a f m-demodutator 
having a considerably reduced number of noise peaks in the demodulated signal and therefore a better signal to noise 
ratio (SNR) below the fm-threshold than the state of the art fm-demodulators. 
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[0021 ] Figure 3 shows an enlarged view of the unwrapped phase <Pu(kT) with a carrier to noise ratio of 1 0dB of a fre- 
quency rrodulated stereo-multiplex signal, this time having a rising slope. This unwrapped phase signal <Pu{kT) has a 
phase junrp of '2n at the sample k = 0. 

[0022] Rgure 4 shows a principal flowchart how the unwrapped phase (Pu(kT) is calculated from the modulo 2n divided 

5 phase <p(kT) according to the present invention. After the start In step S1 it is checked in step S2 whether the difference 
of the phase (p(kT) at sample k and the phase <p{{k - 1) T) at the preceding sample to sample k exeeds n. If this is not 
the case, then it will be checked in step S3 if this difference is below -jl If this is also not the case. i. e. if said difference 
is within the range [-ii. .... n], then the unwrapped phase (f^{kT) is set to the phase <p(kJ). Then the algorithm is set forth 
with step S2 for the next sample k = k + 1 . 

10 [0023] If. on the other hand, it is determined in step S2 that said difference inbetween the phase 9(kT) at sanrple k 
and the phase 9((k - 1 ) T) of the preceding sample of sample k exeeds n then 2n will be subtracted from the actual sam- 
ple of the phase <p(kT) and all following samples of the phase in a step S4 before the algorithm is set forth with step S6. 
[0024] In the same manner 2n will be added to the actual sample of the phase <p(kT) and all following samples of the 
phase in a step S5 if the difference in-between the actual sample of the phase (p(kT) and the preceeding sample of the 

15 phase <p((k - 1) T) falls below -n in step S3 t>efore the algorithm will be set forth with step S6. 

[0025] The phase jumps of the unwrapped phase <Pu(kT) are detected and eliminated with the correction unit 5b 
shown in figure 1 which serves as detection and elimination block that checks the unwrapped phase signal <Pu(kT) for 
positive phase jumps at the falling slope and negative phase jumps at the rising slope in combination wrth an interpola- 
tion before and after the possible phase jump. If a phase jump is delected, the block eliminates the phase jump by add- 

20 ing or subtracting 2n to the following samples of the unwrapped phase signal <Pu(kT). 

[0026] The first criterion for a detection of a phase jump is the rising or tailing slope. This criterion can be determined 
in combination with an interpolation before and after the test sample to be checked for the occurance of a phase jump. 
This interpolation is explained in the following. 

[0027] The interpolation makes use of the spectral distribution of the frequency modulated signal in combination with 
25 the high sampling rate of the unwrapped phase signal <Pu(kT)- In case of fm-stereo-broadcasting. the stereo-multiplex- 
signal s(t) contains mainly low frequency components. The sampling rate of the demodulator has to be chosen high, e. 
g. 1/T > 300 kHz. due to the bandwidth of the baseband signal Or{kT). Therefore, the phase jump has a stronger gradi- 
ent than the absolute phase <Pc(kT) and can be detected. 

[0028] As mentioned before, figure 3 shows the unwrapped phase <Pu(kT) of the fm-demodulaled signal with a phase 
30 jump of '2k at sample k = 0. For the Interpolation at sample k the gradient before k. named grad^ef <Pu(kT). and after the 
sample k, in the following named grad^ft <Pu(kT), has to be calculated. This should be done by averaging the gradient of 
more than one sample, because the unwrapped phase signal <Pu{kT) is disturbed by noise. Also, the samples in the 
direct neighbourhood of the phase jump should not be used for interpolation, because these samples are disturbed by 
the phase jump, as it can be seen in figure 3. Therefore, the gradient before k and after k is calculated by the following 
35 formula: 

^ yu((*^-o)T)-cPu((k-o-p)T) 
gradi^fikl) = (7a) 
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. yu((k-KH r-1)T)-<Pu((k+r-1)T) 
grad^fflkJ) = 



wherein p and q denote the number of samples used for the averaging and o and r denote the minimal distance of the 
45 samples used for calculating the gradient from k. If for example the gradient before and after k is calculated by averaging 
p=q=8 samples beginning with a distance form k of o=r=4 samples, the gradient will be calculated with the following 
formula: 



so gradt^fikl) « 



yu((»<-^)T)-yu(('^-'*^)"0 (7b) 



^ (PuMk^li m-cpuUk-i^)!) 
grad ^ff[kT) g 

[0029] The gradient (7) Is now used for the interpolation. The en-or of the interpolation of the unwrapped phase <Pu(kT) 
is calculated at the sanple k with the gradient before k with the following equation: 
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<P errt>e/, .(kT) = <P u((k-d-1)T) + (d+e+1) • gradt^fikT) - (p u((k+e)T) (8a) 

wherein d and e denote the distance of the interpolated sanples from the actual sample k. In the following, d and e are 
5 supposed to be the same values (d= e ). txjt this is not necessary. 

[0030] Since the sanples in the direct neighbourhood of the phase jump at sample k are disturbed and shouki there- 
fore not be used for the interpolation of the unwrapped phase, d and e should be choosen ^ 2. 
[0031] For improving the reliability of the interpolation error, the interpolation should be done for more than one sam- 
ple. The following example shows the interpolation en'or (Perrbeid,eV<^ 6-e=2 . 3, 4. 5: 

10 

<P^rrt»/^2(kT) = q)u((k-3)T) + 5 • grad ^^^(kT) - (pu((k+2)T) (8b) 
<Perrt>ef33(kT) = <Pu((k-4)'n + 7 • grad^,(kT) - q)u((k+3)T) 

15 

<Peir6e/4.4(k'n = <Pu((k-5)T) + 9 • grad ^,(kT) - (pu((k+4)-0 
<P errbefj^'^ = <P u((k-6)T) + 11* orad ^,(kT) - <p u((k45)T) 

20 

As another point for approving the reliability of the interpolation, It should also be done with the gradient after the phase 
jump. Then, the phase error <Perram,g(kT) will be calculated by the following interpolation: 

(kT) = i>u((k-M)"n + (f+g+1) • gradaff(kT) - <Pu((k+g)T) (9a) 

wherein f and g denote the distance of the interpolated samples from k. f and g are again supposed to be the same 
value (f=g ), even if this is not necessary. 
30 [0032] Because of the unreliability of the neighbour samples of the phase jump at K the minimal distance of the inter- 
polated samples from k should again be choosen f, g ^ 2. The following example shows the interpolation error 
<P«7aftf.g(kT) for f=g= 2 , 3. 4. 5: 

<P«^afr,,(kT) = ^>u{(k-3)T) + 5 • grad^fr(kT)-q)u((k+2)T) (9b) 

35 

9e,ra/f,33(kT) = <Pu((k-4)T) + 7 • grad ^(WV) - cpu((k+3)T) 
<P«r«ft,.,(kT) = <Pu((k-5)T) + 9 • grad^(kT) - qpu((k+4)T) 
^eirafts^ (kT) = 9u((k-6)T) + 11 • grad^,,(kT) - <Pu((k+5)T) 

45 [0033] The addition of all interpolations results in the total interpolation error: 

^ max ' max 

<Perr(kT) = £ ^ orrtefJ^'^ + L <P e/raft, ^(kT) (10a) 
so * 

= £ <P J(k-d-l)T) + (d+efl) 9radt^f{kT) - q)u((k+e)T) 

d=do 

55 ^ max 

+ 2^qp J(k-f-1)T) + (f+g+1) grad^^(}<V) - <Pu((k+g)T). 



40 
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wherein do and fo denote the the nninimal and dn,ax ^ W denote the maximal distance of the interpolated sanples 
from k In the above examples with do=fo=2 and <in^=^max-^ interpolation is done with 

(d max"do+*max **o+ 2)=8 samples, the addition of all interpolations result in tiie total interpolation error: 

<p^^(kT) = 2 • (9u((k-3)T) + (pu((k-4)T) + 9u((k-5)T)+9u((k-6m^ 
- 2 • (<Pu((k+2)T) + <Pu((k+3)T) + qPu((k+^)'n + <P u((k-^5)T)) + 32 • {g rad t^,{kT) + grad^ikl)) 

= 2 • (q)u({k-3)T) + 9u((k-4)T) + (Pu((k-5)T) + (Pu((k-6)T)) (10b) 
-2 • (9u{(k+2)T) + <Pu((k^^)T) + 9u((k^^)"^+9u{(k+5^^^ • (q>u((k-4)T) 
- <Pu((k-12)T) + (Pu((k4.1l)T) - 9u((k+3)'n) 



[0034] In case of an ideal phase jump, the total interpolation error can t>e calculated by the following equation: 

<P e/r^e- = 2ii • (d max<i 0+^ max 0+2) (1 la) 

[0035] In the above example with d o=*o=2 and d max =^ max =5 interpolation error would result in 

20 9^^.^^ = 271 - 8 = 50.27 {11b) 

[0036] On the other hand, practical experiments have proven, that an interpolation error threshold of 

Verr.tfir-I^ • ^Perr^ 02) 



25 



results in good performance at a samplingrate l/T of approximately 410 kHz. For the above exanrple. <Perr,thr = 26 was 
choosen. 

30 [0037] As desaibed before, the unwrapped phase <pu(kT) has to be checked for phase jumps. Since the phase jump 
with the value +2ti appears only at the falling slope of the unwrapped phase signal qpu(kT). the phase jump can be 
detected in the following way. The detection of the falling slope can be done with the gradient of the unwrapped phase. 
In the following, the gradient before kT gradt^^kT) is chosen for the check of this condition (gradt^^kT)<'gradttirY it is 
also possible to chose the gradient after kT (gradsft{kT)<'gradf^r) o"" a combination of botii {gradt^^kJ) + gradgfflkJ) <- 

35 gradthr)* wherein gradfhr is a threshold value preferrabty equal to zero or a small value near to zero. 

((pu(kT) - 9u((k-l)T > <Pthr2(kT)) (13) 
AND (gradbe/(kT) < -gradthr) 

40 

AND (VerrCkT) < -<^err,thr) 



45 [0038] The p)hase jump Is not always such a sharp jump as can be seen from figures 2 and 3. Some phase jumps 
extend over more than one sample so that the following condition is included in the check for the phase jumps according 
to the present invention: 

q>u(kT) - 9u((k-2)T) > <p^,2(kT) (14) 

so 

[0039] With equations (13) and (14) it follows for the condition to detect a phase jump at the falling slope of the 
unwrapped phase signal (Pu(kT): 

55 
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H9u(kT) - 9u((k-l)T) > q>thri(kT)) OR (<pu{kT) - <pu((k-2)T) > q>thr2(kT))l 
AND IgradbeA^'^ < -gradthr) 
AND (<perr(kT) < -q>err.thr) 



(15) 
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[0040] The condition for the rising slope of (Pu(i^T) can be desaibed in a similar way: 



[(<pu(kT) - ipu((k-l)T) < -cpthri(kT)) OR (<pu(kT) - <pu((k-2)T) < -<Pthfi2(kT))l 



(16) 



AND (gradbe/^T) > gradthr) 

AND (<PeiT(kT) > q>err.thr) 



20 [0041 ] The check for the rising slope gradtyefi<T)>gradt^r ca" afla*" be replaced by 9raclsff{kJ)>gradt^r iQf'^^^befi^'^ 
+ 9rad^f^\a))>grad^t,r' 

[0042] Practical experience have proven that the tresholds %>thri and ^thr2 a"cl gradfhr denoted in the above formulas 
(15) and (16) should be chosen at a sampling rate 1/T of approximately 410 kHz as follows: 

25 <p^, «1,5 (17) 



[0043] The flow chart of the detection and elimination is shown in figure 5. This flow chart includes the above-dis- 
cussed check for positive phase jumps at the falling slope of <Pu(kT) and negative phase jumps at the rising slope of 
<Pu(kT) in combination with an interpolation before and after the phase jump. If a phase jump is detected, the block elim- 
inates the phase jump by adding or subtracting 2n to the actual sample and the following samples of the unwrapped 
35 phase signal <Pu(kT). 

[0044] After the start in step SIC, condition (1 6) will be checked in step S1 1 . If this condition is not satisfied, condition 
(15) will be checked in step S12. If this condition is also not satisfied, i. e. if no phase jump has been detected, the cor- 
rected phase <Pe(kT) will be set to the value of the unwrapped phase (Pu(kT) in step SI 5. Thereafter, the algorithm is 
repeated beginning from step S1 1 with the next sample k = k + 1 . 
40 [0045] If, on the other hand, it is in step S1 1 determined that condition (16) is satisfied then 2k will be added to the 
actual sample of the unwrapped phase cpu(kT) and all following samples before step SI 5 is performed. Likewise, if it is 
determined in step 812 that condition (15) is satisfied then 2k will be subtracted from the actual sample of the 
unwrapped phase (Pu(kT) and all following samples before step Si 5 is performed. 

[0046] Figure 6 shows the unwrapped phase cpu(kT) and the phase with eliminated phase jumps <Pe(kT) both with a 
45 carrier to noise ratio of 10 dB after the correction unit 5b. It can clearly be seen that the phase with eliminated phase 
Jumps (Pe(kT) corresponds better to the absolute phase <Pc(kT) with infinite carrier to noise ratio shown in figure 2 then 
the unwrapped phase <^J}CV) without the detection and elimination of the 2n phase jumps. 

[0047] The phase with eliminated phase jumps (PeC^^'O can now be differentiated without the noise caused by the fm- 
threshold. This means that a high gain in the signal to noise ratio of the fm-demodulated signal below the fm-threshold 
so can be achieved. 

[0048] In a second embodiment according to the present invention it is taken into consideration that it is not very easy 
to implement the adding or subtracting of 2n from the actual sample and all following samples of the phase or the 
unwrapped phase, i. e. it is not in an easy way possible e.g. to solve the equation: 

55 fori = kto« (18) 



grad^, « 0 



30 



<p(iT) = <p(iT) - 2x1 
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end 



10 



15 



shown in step S4 o1 figure 4. 

[0049] Therefore, according to the second embodiment of the present invention, the window of calculating the 
unwrapped phase cpu(kT) is limited to the values that are needed for the algorithm. In the akx>ve examples the window 
is limited to the range (k-o-p. ..k-Kj+r-1 ]. e. g. Ik-12...k+1 1 ]. The calculation of the unwrapped phase <Pu(kT) is done with 
the helper variable w, which contains information atxnit the addition or sut>tractton of 2n to the modulo 2n divided phase 

f-l <p(kT) - <pak-l)T) > n (19) 

w(kT) = ] 1 for 9(kT) - 9((k-l)T) < -n 
I O else 



[0050] With this helper variable w(kT) the unwrapped phase (po(kT) can be calculated in a window beginning at sam- 
ple (k-N)T in the following way: 

^ (Pu((k-N}T) = <p((k-N)T) (20a) 

q)u((k.N+l)T) = 9((k-N+l)T) + 2Tcw((k-N+l)T) 

9u((k-N+v)T) = q>((k-N+v)T) + 2:1 1 w((k-N+i)T), 

wherein N denotes the beginning of the window for unwrapping the phase <p(kT) and v denotes the distance of the sam- 
30 pie (k-N+v)T from the beginning (k-N)T of the window for unwrapping the phase. 
[0051] The above example with N = = 12 follows 

<Pu((k- 12)T) = <p((k- 12)T) (20b) 
3^ (Pu((k- 1 1)T) = 9((k- 1 1)T) + 2icw((k- 1 1)T) 

<Pu((k-lO)T) = 9((k-10)T) + 27c(w((k-ll)T) + w((k-10)T)) 

^ <Pu({k-12+v)T) = (p((k-12+v)T) + 2jc X w(k.l2+l). 



45 [0052] With this equation follows for the example of equation (7b): 

grad ^,{kT) = | (9u((k-4)T) - (p^((k-12)T)) (21) 

= 5 (<P((k-4)T) - <p((k-12)T) + 2ii[w((k-1l)T)+w((k-10)T)+...+w{(k-4)T)]) 
so ^ 

[0053] Converting equation (1 Ob) into another representation results in 

55 
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"^^^ = (p J(k-3)T) + (pJ(k-4)-n + 1>u((k-5)T) + q)u((k-6)T) 
- J(k+2)T) + <p^,((k+3)T) + 9u((k+4)T) + ^.^((k+S)!)) + 16 • grad^„(kT) + 16 • gradi,^,m (22) 
= q>u((k-3)T) + *u((k-4)"D + <Pu((k-5)T) + iPu((k-6)T) - (ip^Mk+Z)!) + <Pu((k+3)T) + (Pu((k+4)T) + 
Vuii^+5)J))+2 • (q>u((k+11)T)-9u((k+3)T)) + 16 • grad ^,(kT) 

[0054] With another window definition for the calculation of the unwrapped phase with N = 6 

<Pu((k-6)T) = <p((lt-6)T) (23) 
q>u((k-5)T) = (p((k-5)T) + 2jnv(k-5) 



q>u((k-6+v)T) = <p((k-6+v)T) + 2w Z wlk-6+v) 



follows: 



9err(KT) 



= <p((k-3)T) + i)((k-4)T) + q>((k-5)T) + qp((k-6)-n + 2 • (p((k+11)T) 



2 

- (p((k+2)T) - 3(p((k+3)T) - (p{(k+4)T) - (p((k+5)T) + 2n • {-w((k-5)T) - 2w((k-4)T) - 3w((k-3)-0 - 4w((k-2)T) (24) 
- 4w((k-1)T) - 4w(kT) - 4w((k+1)T) - 4w((k+2)T) - 3w((k+3)T) + w((k+5)T) + 2w((k+6)T) + 2w((k+7)T) 
+ 2w({k+8)T) + 2w{(k+9)T) + 2w((k+10)T) + 2w((k+1 1)T)} + 16 • grad t^f{kT) 

[0055] The equations of the conditions shown in equation (15), (16) and figure 5 can be replaced with the following 
equations: 

9u(kT) - <Pu((k-1)T) = q>{kT) - <p((k-1)T) + 27iw(kT) (25) 

<Pu(kT) - q>u((k-2)T) = <p(kT) - <p((k-2)T) + 27i (w(kT) + w((k-1)T)) 

[0056] The f bw chart of the above example that is shown in figure 5 can be simplified with equations (21). (24) and 
(25) as follows when the condition determined in step S1 1 is replaced by the following condition: 



(26) 



( q)(kT) - <p((k-l)T) + 2icw(kT) < -<pthri OR > 
V q>(kT) - <p((k-2)T) + 2n {w(kT) + w((k-l)T)) < -9uir2 J 
AND gradbef > gradthr AND Cpcrr > <Pcrr»thr 

the condition determined in step SI 2 Is replaced with the following condition 

( (p(kT) - (p{(k-l)T) + 27rw(kn > <pthri OR ^ (27) 

<p(kT) - (p((k-2)T) + 2n (w(kT) -h w((k-l)T)) > cpthr2 J 
AND grad}ief< -gradthr ANJ> cperr < -<Pcrr.thr 

and the differentiation for generating the frequency demodulated signal s(kT) is included in the following equation: 

s(kT) = (p(kT) - <p((k-1)T) + 2iiw(kT) k = k + 1 (28) 
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in step S15. 

[0057] Rgure 7 shows the flow chart corresponding to the (xocess performed in the second emtxxJiment according 
to the present invention. After the beginning of the process in step S20 the condition (26) is determined in step S21 . If 
this condition is not satisfied then condition (27) is checked in step S22. If this condition is also not satisfied then the 
5 differentiation for generating the frequency demodulated signal s(kT) is performed in step S25 according to equation 
(28) before the process is repeated with the next sample k = k + 1 beginning from step S21 . 

[0058] If. on the other hand, it is determined in step S21 that condition (26) is not satisfied then the variable w(kT) will 
be set to w(kT)=w(kT)+ 1 in step S23. Corresponding thereto. w(kT) will be set to w(kT)=w(kT)-1 in step S24 when it 
is determined In step S22 that condition (27) is satisfied before step S25 is performed. 
10 [0059] Figure 8 shows a direct realization of the digital threshold extension unit 5 and the differentiation unit 6 accord- 
ing to the second embodiment of the present invention. This block diagram is realized with delay elements of order T, 
multipliers, adders, subtractors and two comparing units. The block diagram is directly derivable from the above equa- 
tions and conditions (19) to (25). 

[0060] Figure 9 shows the signal to noise ratio of a frequency demodulated music signal with (solid line) and without 
15 (broken line) threshold extension according to the present invention. The simulation was done with parameters as they 
are used for fm-audio broadcasting, so that the frequency deviation is A<o = 2rt • 75 kHz and a deemphase tater is 
included at the demodulator output. The fm-threshold can be seen at approximately CNR = 1 7 dB. The gain in the signal 
to noise ratio caused by the treshold extension is approximately 9 dB. The gain in the SNR is dependent on the fre- 
quency deviation of the modulated signal. The higher the frequency deviation Ao> the higher the gain in SNR. 

20 

Claims 



1 . Method for digitally demodulating a frequency modulated signal (afm(t)). comprising the following steps: 

25 a) calculating a phase (<p) of a digital conr^lex baseband signal Cj{kT) of the received frequency modulated sig- 

nal iaa^)): 

b) calculating a unwrapped phase (<pj of the phase ((p); 

c) detection and elimination of 2n phase jumps in the unwrapped phase (<Pu) to generate a conected phase 
(<Pe): and 

30 d) differentiation of the corrected phase ((pe). 

2. Method according to daim 1, characterized in that said step b) comprises the following steps for every incoming 
sample of the phase (9): 

35 b1) subtracting 2n from the actual sample of the phase ((p(kT)) and a predetermined number (I) of following 

samples of the phase (9((k+1)T)...<p((k+ l)T)) if the difference in-between tiie actual sample of the phase 
(<p(kT)) and the preceding sample of the phase (9((k-l)T)) exceeds ji; and 

b2) adding 27i to the actual sample of the phase ((p(kT)) and a predetermined number (I) of following samples 
of the phase (<p((k+1)T)...(p((k+l)T)) if the difference in-between the actual sample of the phase (q>(kT)) and 
40 the preceding sample of the phase (9((k-1 )T)) fails below - n, 

3. Method according to claim 1 or 2, characterized in that said step c) comprises the following steps for every sample 
of the unwrapped phase (qpj: 

45 c1 ) adding 2/1 to the actual sample of the unwrapped phase ((Pu(kT)) and a predetermined number (1) of follow- 

ing samples of the unwrapped phase (9u((k+1)T)...(Pu((k+l)T)) if the cfifference in-between the actual sample 
of the unwrapped phase (9u(kT)) and tiie preceding sample of the unwrapped phase (9u((k-1 )T)) falls below a 
negative first threshold value { vthri) and the gradient {gradt^f{kT)) of the unwrapped phase before the actual 
sample of the unwrapped phase (<Pu(kT)) or tfie gradient {gradaft(kT)) off the unwrapped phase after the actual 

so sample of the unwrapped phase ((Pu(kT)) or a comtMnation of bofh is higher than a gradient threshold {gradi^) 

: and 

c2) subtracting 2tc from the actual sample of the unwrapped phase (<Pu(kT)) and a predetermined number (I) of 
following samples of the unwrapped phase (<Pu((k+1)T)...<Pu((k*l)T)) the difference in-between the actual 
sample of the unwrapped phase (qPu(kT)) and the preceding sample of the unwrapped phase (<Pu((k-1)T)) 
55 exceeds tiie positive first threshold value (9 thri) and the gradient (gradt^^l)) of the unwrapped phase before 

the actual sample of the unwrapped phase (9u(k"0) or the gradient igrad^ni^T)) of the unwrapped phase after 
the actual sample of the unwrapped phase((Pu(kT)) or a combination off both falls below the negative gradient 
threshold irgrad^). 



11 



BNSOOCID: <EP_<»4005eA1 J_> 



EP0 940 958A1 



m 



4. Method according to anyone of claims 1 to 3. characterized in tliat said step c) comprises the following steps for 
every sample of the unwrapped phase (cpy): 

cl) adding 2n to the actual sample of the unwrapped phase (<Pu(kT)) and a predetermined number (I) of follow- 
5 ing samples of the unwrapped phase (<Pu((k+1)f)...q) y((k+l)T)) If the difference in-between the actual sample 

of the unwrapped phase (9u(kT)) and the preceding sample of the unwrapped phase (<Pu((k-1)T)) falls not 
below a negative first threshold value (-<Pthri) the difference In-between the actual sample of the unwrapped 
phase ((Pu(kT)) ^nd the sample before the preceding sample of the unwrapped phase (<Pu((k-2)T)) falls below a 
negative second threshold value {'<Pthr2)* ^rid the gradient {gradj^fikT)) of the unwrapped phase before the 
10 actual sample of the unwrapped phase (<Pu(kT)) or the gradient (gradani^T)) of the unwrapped phase after the 

actual sample of the unwrapped phase (q>u(kT)) or a combination of both is higher than a gradient threshold 
igradthr) : and 

cll) subtracting 2k from the actual sample of the unwrapped phase (<Pu(kT)) and a predetermined number (I) of 
following samples of the unwrapped phase (q>u((k+1)T)-- <Pu((k+l)T)) if the difference in-between the actual 

15 sample of the unwrapped phase (9u(kT)) and the preceding sample of the unwrapped phase (qpu((k-1 )T)) does 

not exceed the positive first threshold value (^Pthri) l^^t the difference in-between the actual sample of the 
unwrapped phase (<Pu(kT)) and the sample before the preceding sample of the unwrapped phase (<Pu((k-2)T)) 
exceeds the positive second threshold value {<Pthr2)^ and the gradient {gradt^^fikT)) of the unwrapped phase 
before the actual sample of the unwrapped phase (<Pu(kT)) or the gradient {gradsftikJ)) of the unwrapped phase 

20 after the actual sample of the unwrapped phase (<Pu(l^^) or a combination of tK>th falls below the negative gra- 

dient threshold {-grad^f) . 

5. Method according to anyone of claims 2 to 4, characterized in that said predetermined number (I) of following 
samples are all following samples. 

25 

6. Method according to claim 1 . characterized in that said step b) comprises the following steps for every incoming 
sample of the phase (qp): 

b1) setting a helper variable (w(kT)) to -1 if the difference in-between the actual sample of the phase ((T>(kT)) 
30 and the preceding sample of the phase (<p((k-1)T)) exceeds n\ 

b2) setting the helper variable (w(kT)) to 1 if the difference in-between the actual sample of the phase ((p{kT)) 
and the preceding sample of the phase (<p((k-1)T)) falls below - n\ and 

b3) setting the helper variable (w(kT)) to 0 tf the difference in-between the actual sample of the phase ((p(kT)) 
and the preceding sample of the phase ((p((k-1)T)) lies in the interval inbetween -n and ti; 

35 

7. Method according to claim 6. characterized in that said step c) comprises the following steps for every sample of 
the unwrapped phase ((pj: 

cl) adding 1 to said helper variable (w(kT)) if the difference in43etween the actual sample of the unwrapped 
40 phase (q)u(kT)) and the preceding sample of the unwrapped phase (<Pu((k-1)'n) falls below a negative first 

threshold value {-(^thri) and the gradient {gradbef(\CV)) of the unwrapped phase before the actual sample of the 
unwrapped phase (<Pu(kT)) or the gradient {grad^ftiVCY)) of the unwrapped phase after the actual sample of the 
unwrapped phase (<Pu(kT)) or a combination of both is higher than a gradient threshold {gradthr) \ and 
c2) subtracting 1 from said helper variable (w(kT)) if the difference in-between the actual sample of the 
45 unwrapped phase (<Pu(kT)) and the preceding sample of the unwrapped phase (qPu((k-1)T)) exceeds the posi- 

tive first threshold value {^thri) and the gradient {gradt^i(kT)) of the unwrapped phase before the actual sample 
of the unwrapped phase (<Pu(kT)) or the gradient (grao(a/f(kT)) of the unwrapped phase after the actual sample 
of the unwrapped phase ((Pu(kT)) or a combination of both falls below the negative gradient threshold (-grader)- 

so 8. Method according to claim 6 or 7. characterized in that said step c) comprises the following steps for every sample 
of the unwrapped phase ((pj: 

cl) adding 1 to said helper variable (w(kT)) if the difference in-between the actual sample of the unwrapped 
phase (<Pu(kT)) and the preceding sample of the unwrapped phase (qPu((k-1)T)) falls not below a negative first 
55 threshold value {-<Pthri) but the difference in-between the actual sample of the unwrapped phase (<Pu(kT)) and 

the sample before the preceding sample of the unwrapped phase (<Pu((k-2)T)) falls below a negative second 
threshold value (-(p^,,^). and the gradient (drad^/(kT)) of the unwrapped phase before the actual sample of the 
unwrapped phase (q>u(kT)) or the gradient (grada^kT)) of the unwrapped phase after the actual sample of the 
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unwrapped phase ((PuCkT)) or a combination of both is higher than a gradient threshold {grad^) : and 
cll) subtracting 1 from said hefper variable (w{kT)) if the difference in-l>etween the actual sample of the 
unwrapped phase (<Pu(kT)) and the preceding sample of the unwrapped phase (q)u((k-l)T)) does not exceed 
the positive first threshold value {^mti) but the difference in-l>etween the actual sample of the unwrapped 
phase ((Pu(kT)) and the sample before the preceding sample of the unwrapped phase (qPu((k-2)T)) exceeds the 
positive second threshold value (tp^,^). and the gradient (gradbe^kJ)) of the unwrapped phase before the 
actual sanple of the unwrapped phase <Pu(kT)) or the gradient (gradsft(VJ)) of the unwrapped phase after the 
actual sample of the unwrapped phase ((PuCkT)) or a combination of both falls below the negative gradient 
threshold {-gradff^) . 

Method according to anyone of claims 6 to 8. characterized in that the unwrapped phase (<Pu(kT)) is calculated in 
a window beginning at sample (k-N)T in the following way: 

9u((k-N)T) = (p((k-N)T) 

<Pu((k-N+l)T) = 9((k-N+l)T) + 2iiw((k-N+l)T) 



q>u((k-N+v)T) = (p((k-N+v)T) + 27i I w((k-N+i)T), 

20 



wherein N denotes the beginning of the wirxlow for unwrapping the phase (q>(kT)) arxJ v denotes the distance of the 
.sample (k-N+v)T from the beginning (k-rsf)T of the window for unwrapping the phase. 

25 

10- Method according to anyone of claims 3 to 5 or 7 to 9. characterized in that said step c1) or/and cl) is only per- 
formed if the phase error ((Pefr(kT)) of the unwrapped phase at the actual sample of the unwrapped phase ((Pu(kT)) 
exceeds a positive phase error threshold value {<Perr,thr) and said step c2) or/and cll) is only performed if the phase 
error (<Pe„^kT)) of the unwrapped phase at the actual sample of the unwrapped phase (<Pu(kT)) fells below the neg- 
30 ative phase error threshold value (-(Pen^tfir)- 

1 1 ■ Method according to anyone of claims 3 to 5 or 7 to 1 0, characterized in that said first threshold value {<Pthri)' said 
second threshold value ((Ptht^) and said gradiend threshold (gradt^r) are determined by experiments depending on 
the sampling rate (1/T). 

35 

12. Method according to daim 1 1 , characterized in that said first threshold value ((Pthresi) equals to 1 .5. said second 
threshold value ((Pthres2) equals to 2 and said gradiend threshold {gradthr) equals to 0 when said sampling rate (1/7) 
is in the range of 410 kHz. 

40 13. Method according to anyone of claims 3 to 5 or 7 to 12. characterized in that said phase error (<Perr(kT)) of the 
unwrapped phase at the actual sample of the unwrapped phase ((Pu(kT)) is calculated based on the following equa- 
tion: 

^max ' max 

^ <Perr(kT) = £ <Pent^iJ^'^ + L 9erraft.g{kT) 

d«do Mo 

= ^ 9u((k"d-l)T) + (d+e+1) grad ^,{kT) - (Pu((k+e)T) 

50 , 

' mat 

+ L <Pu({k-M)T) + (f+9^-1) grad^^lkl) - 9u({k+g)T). 

1=1 0 



55 wherein grad^^fikT) corresponds to the gradient of the unwrapped phase before the actual sample of the 
unwrapped phase (qpJkT)) and grada^kl) corresponds to the gradient of the unwrapped phase after the actual 
sample of the unwrapped phase (<Pu(kT)). do and fo denote the minimal and dmax and f ^ax denote the maximal dis- 
tance of the interpolated samples from k. and d. e. f and g denote the distance of the sanples used for the interpo- 
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l^vort from k. 



14. Method according to claim 13, characterized in that do = fo = 2,d^ax = *max = S e = d . g = f . 

5 15. Method according to anyone of claims 3 to 5 or 7 to 14, characterized in that said gradient {gradbeA^V) the 
unwrapped phase before the actual sample of the unwrapped phase ((Pu(kT)) is calculated by averaging p samples 
of the unwrapped phase ((pj without those in the direct neighbourhood of the actual sample of the unwrapped 
phase (<Pu(kT)) according to the following equation, the minimal distance to the actual sample k is o, and said gra- 
dient {gradsftikT)) of the unwrapped phase after the actual sample of the unwrapped phase (<Pu(^'0} is calculated 

10 by averaging q samples of the unwrapped phase (cpu) without those in the direct neighbourhood of the actual sam- 
ple of the unwrapped phase ((Pu(kT)) according to the following equation, the minimal distance to the actual sample 
kis r: 

15 9^^d ijgfikT) = 



20 



grad ^^{kT) = 



P 

9u({k4q+r-1)T)-^J(k+r-1)T) 



q 

16. Method according to claim 15, characterized in that o = r = 4 and p = q = 8 . 



17. Method according to anyone of claims 10 to 16, characterized in that said phase error threshold value {<Perr,thr) 
determined by experiments depending on the sampling rate (l/T) and a good start value for this experiments is: 

25 

18. Method according to claim 1 7, characterized in that phase en-or threshold value {<PQrr,thr) equals to 26 when said 
30 sampling rate (1/T) is in the range of 410 kHz. 

19. Method according to anyone of claims 1 to 18. characterized in that said calculation of the phase of the digital 
complex baseband signal CjCkT] is performed using the CORDIC algorithm. 

35 20. Device for digitally demodulating a frequency modulated signal (afni(t)), comprising: 

a) a transformation unit (4) to calculate a pliase {<p) of a digital complex baseband signal cj{kT) of the received 
frequency modulated signal (afm(t)}; 

b) a phase unwrap unit (5a) to calculate an unwrapped phase (cpj of the phase (<p) output of the transformation 
40 unit (4); 

c) a correction unit (5b) to detect and eliminate 2n phase jumps in the unwrapped phase(9u) to output a 
corrected phase (9©) ; and 

d) a differentiatk>n unit (6) to differentiate the corrected phiase (q>e). 

45 21. Device according to claim 20, characterized in that it is working according to the method defined in one of the 
claims 1 to 19 above. 



so 



55 
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Figuie 4: 
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□ BLACK BORDERS 

□ IMAGE CUT OFF AT TOP, BOTTOM OR SffiES 

□ FADED TEXT OR DRAWING 

□ BLURRED OR ILLEGIBLE TEXT OR DRAWING 

□ SKEWED/SLANTED IMAGES 

□ COLOR OR BLACK AND WHITE PHOTOGRAPHS 

□ GRAY SCALE DOCUMENTS 

□ LINES OR MARKS ON ORIGINAL DOCUMENT 

□ REFERENCE(S) OR EXHIBIT(S) SUBMITTED ARE POOR QUALITY 

□ OTHER: 

IMAGES ARE BEST AVAILABLE COPY. 
As rescanning these documents will not correct the image 
problems checked, please do not report these problems to 
the IFW Image Problem Mailbox. 
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